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Abstract 

In this note, we explain how the algorithms in [11], [5] and [6] im- 
mediately give formulas which can be used for the efficient symbolic 
computation of series expansions to solutions of nonlinear systems of 
ordinary differential equations. As a by product of this analysis, we 
derive formulas relating trees to the coefficients of the series expan- 
sions, similar to the work by Leroux and Viennot [15], [16], [17] , and 
Lamnabhi, Leroux and Viennot [18]. 

1 Introduction 

In this note, we explain how the algorithms in [11], [5] and [6] immediately 
give formulas which can be used for the efficient symbolic computation of 
series expansions to solutions of nonlinear systems of ordinary differential 
equations. As a by product of this analysis, we derive formulas relating trees 
to the coefficients of the series expansions, similar to the work by Leroux 
and Viennot [15], [16], [17], and Lamnabhi, Leroux and Viennot [18]. In this 
section and the next, we follow the exposition in [14]. This is an extended 
abstract: a complete version will be published elsewhere. 

We now describe the basic idea of how trees can be used to organize 
computations involving vector fields following [8] and [7]. For background 
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material on Hopf algebras, see [19]. Consider a control system 

i(t) = E l (x(t)) + u 1 (t)E 2 (x(t)) + u 2 (t)E 3 (x(t)), x(0) = x°eK N , (1) 

where E \ , E 2 and E 3 are vector fields defined in a neighborhood of x° £ R^ 
and t 1 — ► u{(t) are controls. Our goal is to describe a class of algorithms for 
the effective symbolic computation of expressions built from the vector fields 
that describe the local behavior of the control system. These expressions 
include iterated Lie brackets and the generating series of the system. 

The starting point is to assign trees to vector fields as illustrated in Fig- 
ure 1, and then to impose a multiplication on the trees which is compatible 
with the composition of vector fields. Assume that the vector fields Ej have 
the form: 

N 

A = X>?A.. i = 1 . 2 < 3 > j = !,■■■, M, (2) 

1 

where cij are smooth functions on R^ and = d/dx^. Now 
E 2 • E\ = ^2 bj (Djdi) D{ + ^2 bjdiDjDi 


and E 3 • E 2 • E\ is equal to 

E A°l* ) A +£ a 3* (At A< )A+£ a? a?(Dja ? ) D k A 


+Y,a3 k a2< At A A +2 a 3 * a 2 ' ’ ( A O A A+5>£* af ( A ) A A • 

( 3 ) 


Here the sum is for i, j, k = 1, . . . , N and hence involves 0{N 3 ) differentia- 
tions. It is convenient to keep track of the terms that arise in this way using 
labeled trees: we indicate in Figure 2 the trees that are associated with the 
six sums in this expression. 

An iterated Lie bracket such as 


[£ 3 , [E 2 , Ex]] = E 3 E 2 E X - E 3 E X E 2 - E 2 E x E 3 + E X E 2 E 3 (4) 

gives rise in this fashion to 24 trees corresponding to the 24N 3 differenti- 
ations that a naive computation of this expression requires. On the other 
hand, 18 of the trees cancel, saving us from computing 18N 3 terms. We 
are left with QN 3 terms of the form (junk)D^. A careful examination of 
this correspondence between labeled trees and expressions involving the Ej ’ s 
shows that the composition of the vector fields Ej’ s, viewed as first order 
differential operators, corresponds to a multiplication on trees. This multi- 
plication is illustrated in Figure 3. It turns out that this construction yields 
an algebra, which we call the algebra of Cayley trees . 
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2 The algebra of Cayley trees 

In this section, we follow [13] and define a bialgebra structure on spaces of 
trees. The relation between trees and differential operators goes back at 
least as far as Cayley [3] and [4]. Of this literature, the work most closely 
related to the view point taken here is Butcher’s use of trees to analyze 
Runge-Kutta algorithms [1] and [2], 

Let k will denote a field of characteristic 0. By a tree we will mean a 
finite rooted tree. Let T be the set of finite rooted trees, and let be 

the ifc-vector space which has T as a basis. 

We first define the coalgebra structure on k{T}. If t £ T is a tree whose 
root has children s x , . . . , s r , the coproduct A (t) is the sum of the T terms 
t\ ® t 2) where the children of the root of t\ and the children of the root 
of t 2 range over all 2 r possible partitions of the children of the root of t 
into two subsets. The augmentation e which sends the trivial tree to 1 and 
every other tree to 0 is a counit for this coproduct. It is immediate that 
comultiplication is cocommutative. 

We next define the algebra structure on k{T}. Suppose that t\ } t 2 € T 
are trees. Let s\, . . ., s r be the children of the root of t\. If *2 has n + 1 
nodes (counting the root), there are (n + l) r ways to attach the r subtrees of 
t ! which have , . . . , s r as roots to the tree 1 2 by making each s,* the child of 
some node of t 2 . The product t\t 2 is defined to be the sum of these (n + 1) 
trees. It can be shown that this product is associative, and that the trivial 
tree consisting only of the root is a right and left unit for this product. It 
can also be shown that the maps defining the coalgebra structure are algebra 
homomorphisms, so that k{T } is a bialgebra. For details, see [9]. 

The bialgebra k{T } is graded: k{T} n has as basis all trees with n + 1 
nodes. Because the bialgebra k{T} is graded connected, it is a Hopf algebra. 
We summarize the above discussion in the following theorem. 

Theorem 2.1 The vector space k{T} with basis the set of finite rooted trees 
is a cocommutative graded connected Hopf algebra. 

Assume now that each node of the tree (except for the root) is labeled 
with a symbol from the set {E\ ) . . ., Sm}. We can define the product and 
coproduct as above, and, once again, the resulting space is a bialgebra. See 
[9] for details. Let k{CT } denote this algebra. 

Let R denote a subring of the commutative ring of smooth functions on 
R^. We now define an action of the algebra of Cayley trees 

B = k{CT } 
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on the ring /2, making R a B-module algebra, which captures the action 
of trees as higher derivations. This requires that we interpret the formal 
symbols Ej as derivations of R using Equations 2. The action is defined 
using the map 

t J> : k{CT} — ► End/. R y 

as follows: 

1. Given a labeled, ordered tree t with m + 1 nodes, assign the root 
the number 0 and assign the remaining nodes the numbers 1, . . . , m. 
We identify the node with the number assigned to it. To the node k 
associate the summation index p*. Denote (/i i, . . . , jim) by /j. 

2. For the labeled tree t y let it be a node of t } labeled with E^ k if k > 0, 
and let , /' be the children of k . Define 

R{k]fi) = if k > 0 is not the root; 

= , if k = 0 is the root. 

Note that if k > 0, then R(k]fi) (= R. 

3. Define 

N 

4. Extend to all of k{CT} by linearity. 

It is straightforward to check that this action of B on R makes R into a 
5-module algebra. 

We summarize with the following theorem. 

Theorem 2.2 Let R the algebra of smooth functions on R^. Let B denote 
the algebra of Cayley trees k{CT}. Then R is a B -module algebra with 
respect to the action defined by t/>. 

3 Taylor flows for vector fields with polynomial coefficients 

The standard action of the algebra A of differential operators generated by 
Ei, . . ., Em on the algebra of smooth functions R gives R the structure of a 
j4-module algebra. It is easy to relate this /f-module algebra structures on 
R to the ff-module algebra structure defined in the section above. 
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Let 


B 


<f-.A — * 

denote the map sending the generator Ej of the algebra A to the tree con- 
sisting of two nodes: the root and a single child labeled Ej. This map is 
illustrated in Figure 1. Extend <j> to be an algebra homomorphism. Let \ 
denote the map 

A — > Endjt R 

defined by using the substitution (2) and simplifying to obtain an endomor- 
phim of R. We have the following diagram: 

A -* B 

\ 1 ( 5 ) 

Endjt R 

Theorem 3.1 (i) The maps \> <f> an d "P are related by x = *P ° 4>' (**) 
a function f G R and a differential operator p G A. Then 

p - f = <Kp) • /• 

Here the action on the left views R as an A-module algebra, while the action 
on the right views R as B-module algebra. 

The first assertion is proved in [12] and the second assertion follows from 
the first assertion and the definitions. 

To simplify the notation, in the remainder of this section we restrict 
attention to a single derivation F . We will need the following Theorem. 

Theorem 3.2 Assume f 6 R and F G Der (R). Making use of the actions 
defined in Theorem 3.1 , we have: 

(i) For k > 0, 

(F k f)(x) = ^ f(exp(tF)x ) | t =o • 

(ii) If f is analytic near x, then for sufficiently small t, 

OO j_k 

f(exp{tF)x) = ^2f{x-,F k )- 
k=o 

where f(x\F k ) is defined to be ( F k f)(x ). 

(in) If f is analytic near x, then for sufficiently small t, 

/(exp(fF)i) = exp (t<j>(F)) ■ f\ x . 
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Proof. Assertions (i) and (ii) can be found in [20]. Theorem 3.1 and As- 
sertion (ii) then give Assertion (iii). 

Assume now that the derivation F is of the form 

F= £ U»x»D„ 

Ml = l 


Ml ,M2,M3 = 1 

where all except a finite number of the coefficients £ k are zero. 

We let kT{F,x] denote the vector space with basis the set of finite, 
rooted trees, whose nodes are labeled with the symbols F and x, and let 
k<F> denote the free associative algebra generated by the symbol F . We 
define 

<f> : k<F> — > JtT{F,x} 

on the generator F as indicated in Figure 4, and extend cf> to all of k<F> by 
making it an algebra homomorphism. Figure 5 illustrates the multiplication 
in kT {F, x }. 

We define a map 

iP : fcT{F,x} — » Diff(fl) 

as follows; 

1. Given a tree t £ fcT{F, x] with m + 1 nodes, assign the root the 
number 0 and assign the remaining nodes the numbers 1 , . . . , m. 

2. Consider a node with children k ) . k f labeled with an F and 

children labeled with an x. Define 

R{j] A 1 ) — 

if the node is labeled with a F; 

— ’ 

if the node is labeled with a x; 

— ’ * ‘ 

if the node is the root 


(Note that in the first case, is zero, since the F s are constant.) 
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3. Define 


N 

ip(t) = ^2 R(m-,fi)---R{l;ft)R(0;fi). 

Ml»*--,Mm=l 

4. Extend to all of kT{F,x} by linearity. 

We can now state our main theorem. 

Theorem 3.3 Let the derivation F be of the form 

P=f;t«D„+ £ 6£x”D„, 

Ml =1 Ml.M2=l 

Ml .M2,M3 = 1 

where 6£‘ € fc and at most a finite number are nonzero. Then the ith 

component of the nth term in the formal series expansion of the flow 

x(i) = F(x(t)), *(0) = *°, x(t) G R* 


is ^iwen 6y 

0C U ) ' *’1*=*°- 

U 

u;/iere the sum is over all labeled trees u £ fcT{F, x) occuring in the product 
(</>(F)) n . Here the action of a tree u on the function x* is determined by 
viewing R as a B-module algebra using the action just defined above f and 
0(F) is defined as in Figure f 

This follows from Theorems 3.1 and 3.2. 
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